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ABSTRACT 



Aims. The radiance and Doppler-shift distributions across the solar network provide observational constraints of two-dimensional 
modeling of transition-region emission and flows in coronal funnels. These distributions have not, however, been studied in detail and 
we attempt an investigation for a quiet Sun region. 

Methods. Two different methods, dispersion plots and average-profile studies, were applied to investigate these distributions for three 
EUV lines. In the dispersion plots, we divided the entire quiet Sun region scanned by SUMER into a bright and a dark part according 
to an image of Fe xn taken by EIT during the scanning; we plotted intensities and Doppler shifts in each bin as determined according 
to a filtered intensity of Si n. We also studied the difference in height variations of the magnetic field as extrapolated from the MDI 
magnetogram, in and outside network, in the two parts. For the average-profile study, we selected 74 individual cases and derived 
the average profiles of intensities and Doppler shifts across the network. Cases with large values of blue shift of Ne vm were further 
studied. 

Results. The dispersion plots reveal that the intensities of Si n and C iv increase from network boundary to network center in both 
the bright and dark parts. However, the intensity of Ne vm shows different trends, namely increasing in the bright part and decreasing 
in the dark part. In both parts, the Doppler shift of C iv increases steadily from internetwork to network center. The height variations 
in the magnetic field imply a more homogeneous magnetic structure at greater heights and clearly reflect the different magnetic 
structures in the two regions. The average-profile study reveals that the intensities of the three lines all decline from the network 
center to internetwork region. The binned intensities of Si n and Ne vm have a good correlation. We also find that the large blue shift 
of Ne vm does not coincide with large red shift of C iv. 

Conclusions. Our results suggest that the network structure is still prominent at the layer where Ne vm is formed in the quiet Sun, and 
that the magnetic structures expand more strongly in the dark part than in the bright part of this quiet Sun region. Our results might 
also hint for a scenario of magnetic reconnection between open funnels and side loops. 

Key words. Sun: corona-Sun: transition region-Sun: UV radiation-Sun: magnetic fields 



1. Introduction 

Our knowledge of the solar atmosphere, in particular the transi- 
tion region where temperature increases and density drops dra- 
matically with height, depends largely on spectroscopic observa- 
tions. Ultraviolet emission lines can provide ample information 
about the physics of the upper chromosphere, transition region, 
and lower corona. Since the plasma is optically thin for most 
emission lines, we can extract information about the physical 
conditions prevailin g in their sourc e reg ions from t he line pro- 
files (see reviews by iMariskald 1992b and lXial J2003h . 

Before the launch of SOHO (the Solar and Heliospheric 
Observatory), our understanding of the transition region were 
largely based on observations made by early spectrographs such 
as the Naval Research Laboratory (NRL) S082-B EUV spec - 
trograph onboard the Skylab space station (Bartoe e t alj|1977h . 
and NRL High-Resolution Telescope Spectrograph (HRTS) 
flown on some rockets and Spacelab2 dBrueckner et "aLI 1 1977 



iBrueckner and Bartoell 19831; iBrueckner et al.ll 1 986b - These early 
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observations provided much valuable infor mation about the so- 
lar transition region, which was reviewed by Ma riskal dl992b . 

Due to the high spectral and spatial resolution of the SUMER 
(Sola r Ultraviolet Measurements of Emitted R adiation) instru- 
ment dWilhelm et aLlfT^lLemaire et al.fl997b onboard SOHO 
(Solar and Heliospheric Observatory), many more line profiles 
were obtain ed, identified and used in the studies dCurdt et al.l 
1200 ll 12004b . and the Doppler shifts of these EUV spectral 
lines can be measured with an accuracy of about 1-2 km/s 
dBrekke et alJ | l997t i Hassler et all 1 19991: iPeter and Judge! 1 19991; 
Wilhelm et al. I200a iPopescu et alJ 12004b . Besides, another 
spectrograph CDS (Coronal Diagnostic Spectrometer) onboard 
SOHO also has the capability of simultaneous multi-wavelength 
imaging in a numb er of transition region and coronal lines 
dHarrison et alJl995b . Thus, our understanding of the solar tran- 
sition region was improved significantly since 1995 when SOHO 
was launched. 

One of the most prominent features in the chromospher e 
and transition region is the network structure dReevesiri976h . 
which is the upward extension of the supergranular boundary 
at greater heights, and is characterized by clusters of magnetic 
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flux concentrations (Pa tsourakos et al.f i 999). The network man- 
ifests itself as bright lanes in the radiance images of chro- 
mospheric and transition region lines. The network cell has a 
size of 20.000-30.0 00 km dSimon and Leightonlll964l). an d lasts 
for 20 to 50 hou rs dSimon and Leightonlll964tlSchriiverlll997t 
iRaju. et ail 1 19981). The height variation in the network s ize ca n 
be found in iPatsourakos et ail d!999l) . iGontikakis et al.l (120031) . 
an dlTian et al.l d2008bl) . ~ 

iGabriell (I 1976b proposed the first magnetic network model, 
in which the transition region emission originates from funnels 
diverging with height from the underlying supergranular bound- 
ary. However, this model failed to r eproduce the e mission mea- 
sure a t temperatures below 10 5 K (lAthay|[l982l) . iDowdv et alj 
(119861) discovered many fine-scale structures of mixed polarities 
in the photospheric magnetic network; they proposed a modified 
model in which only a fraction of the network flux opens into 
the corona as a funnel shape, while the rest of the network is 
occupied by a population of low-lying loops with lengths less 
than 10 4 km. In this model, the s mall loops are heated internally. 
Based on SUMER observations, iPeterl (1200 ll) suggested a new 
sketch of the structure of the solar atmosphere, which includes 
different structures that dominate at different heights. The most 
intriguing feature in this structure is that there are two types of 
coronal funnels: funnels connected to the solar wind and fun- 
nels that are feet of large loops. It might also be possible that 
the fundamental structure of the solar tra nsition region has not 
yet been resolved by current observations ( Feldmanll983l 119871 : 



iFeldman and Laming||l994] : lDoschek et al 



2004 . 



Another prominent characteristic of the solar transition re- 
gion is the observed redshi ft of transition region lines (see 
a review in iMariskal dl992l) ). The redshift can reach up to 
15 km/s in the network. In a statistical sense, these redshifts 
imply the presence of plasma flows or wave motions in the 
quiet Sun wi th amplitudes that are substantial fractions of the 
sound speed dPeter and Judgell999l) . Different mechanisms such 
as solar wind outflows, spicules, siphon flow through loops, 
nano-flares, explosive events, and downward propagating MHD 
waves have been suggested to explain these pr ominent Dopple r 
shifts of tran s ition r e gion lines (see reviews in IMariskal d 19921) . 
iBrekke et al.l d 19971) . iPeter and Judge! dl999l) . and lXial 420031) 1 
The redshift is most significant in the middle transition region 
and decreases to become a blue shift in the upper transition 
region. The redshift to blueshift transition occurs at electron 
temperatures of about 5 x 10 5 K and the Doppler shift varia - 
tion with temperatu re ca n be found in IPeter and Judge! d!999l) . 
iTeriaca et al.l dl999t) . and lXia et al.l d2004l) . 

There is an increasing consensus that the observed large blue 
shifts in the upper transition region and lower corona are associ- 
ated with solar wind outflows in coronal holes. Relative to past 
studies, a far more accurate rest wavelength of Ne vm (A11Q, 
2s 2 S \/2-2p 2 P3/2) which originates in t he upper transitio n regio n 
and lower corona, were measured by iDammasch et al.l (1 19991) . 
The Doppler shift of this line can be therefore determined well 
and has been intensively studied. Sizable areas of large blue shift 
were found in the Ne vm dopplergrams in coronal holes and con- 
sidered to be a s i gnature of solar wind outflow (lHassler et alj 
1 1 999t iPeterll 1 999t Istucki et al.ll2000t Iwilhelm et al.ll2000l) ; thev 
tended to be la r ger in the darker r egions of coronal holes 
dXia et alj 120031) . lHassler et al.l d 19991) studied the relationship 
between Doppler shift and chromospheric network, and found 
that a larger blue shift was c losely associated with the underly- 
ing chromospheric network. iTu et alj (12005 al) inferred that the 
patches of large Ne vm blue shift were associated with coronal 
funnels reconstructed by a force-free field model. 



In the quiet Sun, large Ne vm blue shifts were also found in 
the ne twork lanes, and co nsidered to be a signature of the solar 
wind (lHassler et al.l 1999|) or just mass supply to quiet coronal 
loops Me et al.ll2007b iTian et alj|2008al) . 

However, the distributions of intensities and Doppler shifts 
of EUV lines across the network, which are required for 2-D 
modeling of transition region emissions and flows in corona l 
funnels, have not been investigated well. lAiouaz et al. 1 d2005bl) 
developed a technique of three-step filter, which was applied to 
an image of Lymann continuum intensity to reveal the chromo- 
spheric network. They binned the filtered continuum intensity 
and calculated the median value of the Doppler shift of Ne vm 
in each intensity bin to produce the so-called dispersion plot. 
They found that the maximum outflow (blue shift) does not ap- 
pear in the center o f the network but rather near network bound- 
ary. lAiouaz! d2008l) extended this work by studying the redshifts 
of transition region lines and found the locations of the maxi- 
mum inflows (redshift) of these lines are network cen ters in the 
quiet S un and network boundaries in the coronal hole. IXia et al] 
(|2004|) studied the correlation between radiance and Doppler 
shift of Ne vm and found that a slightly positive correlation ex- 
ists. They also found that the correlation is not clearly present 
for the C iv line. 

Magnetic field extrapolation is a common way for the solar 
community to st udy the magnetic coupling o f differ e nt solar pro- 
cesses (see e.g. IWiegelmann and Neukirchl d2002l) ; iRuan et ail 
(120081) ). With the assumption that the magnetic field may be con- 
sidered as almos t force-free at heig h ts of about 400 km above 
the photosphere dMetcalf et al.lll995l) . lMarsch et ai] d2004 com- 
bined the extrapola ted magnetic fie ld based on the force-free 
model proposed by Se ehaferl d 19781) with EUV observations to 
study the plasma flows in ARs and detected strong adjacent 
up and down flows associat ed with sunspots. A simil a r method 
was applied successfully by IWiegelmann et alj (12005b . ITu et al.l 
(l2005al) . lTu et all d2005cT), andlMarsch et all d2006l) fo r a coronal 
hole stu dy, in lTu et al.l (l2005bl) . iHe et all d2007l)rand ITian et al.l 
(l2008al) for a quiet sun study, and in ITian et alj d2007b for the 
study of coronal bright points. 

In this paper, we use two different methods, namely the 
method of dispersion plot and average-profile study, to investi- 
gate the distributions of radiance and Doppler shift of three EUV 
lines across the network in a quiet Sun region. 



2. Observation and data reduction 

On 22 September 1996, a middle-latitude quiet-Sun region was 
scanned by slit 2 (1" x 300") of the SUMER instrument from 
00:40 to 08:15 UTC. The raster increment was about 3". Three 
strong EUV lines, Si n (153.3 nm), C iv (154.8 nm), and 
Ne vm (77.0 nm) were included in the selected spectral window 



This d at a set was inten s ively st udied before b y 



(119991) lHassler etalJ dl9 99h IPeterl d2000t). 



Dammasch et al 



Gontikakis et al 



(120031) . rTuet alj d2005bl) . and ITian et all d2008al) . More details 
about this observation can be found in these papers. The selected 
area in our study is 442" x 259" in size, which is the same as 
that considered bv lTian et al.l (l2008al) . 

The standard SUMER procedures for correcting and cali- 
brating the data were applied, including local-gain correction, 
flat-field correction, geometrical-distortion, and dead-time cor- 
rection. 

Each observed spectrum was fitted by a single Gaussian 
curve, complemented by a constant and a linear term describing 
the background. By integrating over a fitted profile, we derived 
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Fig. 1. Upper: The intensity image of Si n, the solid and dashed 
contours represent strong magnetic fluxes with different polari- 
ties (larger than 15 Gauss). Lower: Image of EIT-Fe xn observed 
at 02:30 on 22 September 1996, the red lines mark the bound- 
ary between the bright part and the dark part. The intensities are 
shown in a logarithmic scale. The green and yellow circles mark 
positions of the data points in bins representing the network cen- 
ter in Figs.[3]and |H respectively. 



the total count rate for the profile and thus obtained its inten- 
sity that was then used to construct the intensity image of each 
line. The Doppler shi ft of each line wa s determined by using the 
method described bv lTu et all d2005bl) . 

We used th e magnetogram take n by the Michelson Doppler 
Imager (MDI) (IScherrer et al.l 19951) onboard SOHO on the same 
day at 01:39 UTC. The pixel size of the magnetogram is about 
2". The correction of the magnetogram and the coalignment of 
the magnetogram with the SUME R images were ca rried out by 
applying the method described by Tu et al. (2005b). The linear 
Pearson correlation coefficient between the coaligned magne- 
togram and Si n intensity image is 0.15 (the number of data 
points is 29436 and the corresponding critical correlation co- 
efficient is 0.01, so the magnetic field and chromospheric emis- 
sion correlate very well). The Ne vm intensity image was used to 
coalign the SUMER and EIT images. The uncertainties in both 
coalignments are less than 2". 

We normalized the intensities with respect to the median 
value over the entire area for each line. The intensity image of 
Si n is presented in the upper panel of Fig. [1] where the solid and 
dashed contours represent strong magnetic fluxes of different po- 
larities (higher than 15 Gauss). The network pattern is clear in 



this chromospheric line and almost all the strong magnetic fluxes 
are located in networks. The lower panel of the figure shows 
the sub-reg ion of a full-disk EIT(Extr eme-Ultraviolet Imaging 
Telescope) dDelaboudiniere et al.| [l995) image corresponding to 
the scan area of SUMER. The EIT image was observed in Fe xn 
at 02:30 on 22 September 1996. 

3. The filtered intensity image of Si n 

If the network pattern has a regular shape and the radiance de- 
creases from the center of the network to the internetwork, then 
the radiance can be used to indicate the distance from the net- 
work center. A strong, median, and weak radiance should be 
found close to the network center, network boundary, and the 
internetwork region, respectively. In Fig. Q] we can discern that 
the network structure is very clear and thus the intensity of Si n 
can be used to represent the distance from the network center. 

However, the intensity contrast between the network and in- 
ternetwork for different supergranular cell is different. To re- 
move the local brightenings and reduce th e contrast, we applied 
a three-step filter technique developed by lAiouaz et al. 1 (I2005bl) 
to the intensity image of Si rx. The technique was a pplied instead 
to an in tensity image of Lymann continuum by lAiouaz et al.l 
(I2005bl) . The method can be described in the following way: 
First we applied a median filter to the image. The width of 
the median filter was chosen to be 30", corresponding to the 
mean size of a supergranular cell. The original intensity image 
was then divided by the median-filtered image. Finally, the re- 
sulting image was smoothed over 10", which is the usually ac- 
cepted value of the netwo rk size at transition region temperatures 
dPatsourakos et al.ll 19991) . By using this method, local brighten- 
ings and intensity contrasts of network/internetwork for differ- 
ent supergranular cells were reduced although remain in the fil- 
tered image. We therefore applied two additional iterations of 
this method: the method was reapplied to the image produced by 
the previous filtering iteration. The final filtered intensity image 
of Si ii is presented in Fig. [7] 

In this filtered intensity image of Si n, local brightenings 
have been significantly reduced and the intensity contrasts be- 
tween the network and internetwork for different supergranular 
cells are almost comparable. This filtered intensity was used to 
indicate the distance from network center. A larger, median, and 
smaller value of the filtered intensity of Si n corresponds to a 
location in the network center, network boundary, and internet- 
work region, respectively. Then we are able to complete a statis- 
tic study on the trend of a specific parameter from network center 
to internetwork region, based on this conversion. 



4. Dispersion plots 

Our first method of studying the distributions of radiance and 
Doppler shift across the net work is the method o f di spersion 
plot pr eviously employed by lAiouaz et al.l (l2005bl) and lAiouaz! 
(120081) . We binned this filtered intensity image of Si n and calcu- 
lated the average values of the original intensities and Doppler 
shifts of the included lines in each bin to produce the dispersion 
plots. 

We divided the entire region into two parts of the same area 
according to the intensity of EIT-Fe xn. The two parts are la- 
beled as bright part and dark part. The red lines in the lower 
panel of Fig. Q] indicat e the boundary between the bright and 
dark part. The Fig. 3 of iHe et all d2007l) illustrated spatial rela- 
tions between intersection points of the extrapolated open field 
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lines at 4 Mm, which implies that the open field lines are located 
in the dark part, if we compare it with Fig. Q]of this paper. While 
the bright part is o ccupied by a pop ulation of magnetic loops 
of different heights dTian et al.ll2008al) . Although the correspon- 
dence between the open/closed field lines and the dark/bright 
parts is not exact, a comparison between the statistical study of 
the two parts should highlight a different property of the two 
parts. Another reason why we divided the region into two parts 
is that the network emission of Si n in the bright part is systemat- 
ically higher than that in the dark part. It is therefore more accu- 
rate to present dispersion plots for each part separately to reduce 
the effects of mixing bright internetworks and dark networks. 

To study the relationship between the filtered intensity of 
Si ii (an indication of the distance from network center) and 
the original radiance as well as Doppler shift, we first binned 
the filtered intensity of Si n. The number of data points in 
each bin is shown in Fig. [2] The horizontal axis is given by 
log(Ifut.si nl < Ifiit.si ii >), where If UtSi u is the filtered inten- 
sity of Si n and < 1/ut.si n > is the median value of the filtered 
intensity of Si n over the entire area. It appears that the filtered 
Si ii intensity has a log-normal distribution in both the bright and 
dark part. The sampling steps of the intensity bins were approx- 
imately 0.015. The minimum number of data points per bin was 
chosen to be N„,j„ = 10. If the number of data points in a bin was 
less than 10, the step was changed to 0.030 to make a new bin. If 
the number of data points in the new bin remained less than 10, 
we discarded this bin. 

We then calculated the average values of the original inten- 
sities and Doppler shifts of all three lines and the correspond- 
ing standard deviations in each bin. The dispersion plots for the 
bright and dark part are presented in Figs.[3]and [4] respectively. 
For each plot, the horizontal axis is the filtered Si n intensity 
(log(Ifn,,si nl < Ifiit.si ii >))> which might be considered to rep- 
resent the distance from the network center. The higher the fil- 
tered Si ii intensity, the closer the distance to network center. In 
each of the two parts, we can determine two thresholds for the 
filtered Si ii intensi ty by assuming that the network occupies 1/3 
of the entire area dXia et alj|2004l) and the network center oc- 
cupies 1/25 of the network area; these thresholds are indicated 
by the dashed line representing the outer edge of the network 
boundary, and the dotted line marking the boundary of the net- 
work center in each panel of the figures. In each panel, we can 
consider that the bins to the left of the dashed line are taken from 
the internetwork (cell) region, the bins to the right of the dotted 
line are taken from the network center, and the bins between the 
dashed line and the dotted line represent the network boundary. 

It is unsurprising that the intensities of Si ii and C iv show 
a general increasing trend from network cell to network center, 
since the network pattern is clearly discernible in intensity im- 
ages of chromospheric and transition-region lines. This increas- 
ing trend is obvious within the network, indicating that network 
emission increases steadily from the network boundary to net- 
work center. The result might suggest that the magnetic fluxes 
associated with the EUV emission within the network tend to 
cluster towards the network center. By considering the error 
bars, the chromospheric and transition-region emission appears, 
however, to be distributed uniformly in the internetwork region, 
which suggests that the associated internetwork magnetic fluxes 
are randomly distributed. 

Ne vni emission originates in the upper transition region and 
lower corona. In Figs. [3] and |4j we find that the intensity dis- 
tribution of Ne vni between internetwork and network center is 
similar to those of Si n and C iv in the bright part; in the dark 
part, however, they are not consistent. From the network bound- 



ary to the network center, the intensity of Ne vm shows differ- 
ent trends, increasing in the bright part and decreasing in the 
dark part. Our previous results demonstrated that magnetic struc- 
tures expand throughout the upper transition region and lower 
corona far more sig nificantly with heig ht inside the coronal hole 
than the quiet Sun dTian et alj|2008bl) . And we also found that 
the dark part is associated with open field lines as extrapolated 
from th e photospheric magnetogram dHe et al.ll2.Q07t iTian et ail 
I2008al) . which suggests that the dark part has a similar charac- 
teristic (expanding stronger than the normal quiet Sun) to that of 
the coronal hole. Thus, we suggest the decrease of Ne vm inten- 
sity in the network center is due to a more significant expansion 
of magnetic structures in the dark part. 

The Doppler shift of Ne vm also exhibits a different trends 
between the network boundary and the network center in the 
two parts, if we consider only the average value in each bin. 
In the bright part, the maximum blue shift of Ne vm is close 
to the network boundary, while i t is found in t he netw ork cen- 
ter in the dark part. We note that lAiouaz et al.l d2005bl) claimed 
that the blue shift of Ne vm is strongest at the edge of net- 
work patches in a coronal hole, and the adjacent quiet Sun. They 
conjectured that bipolar fields being swept in from the cell in- 
terior cause enhanced magnetic activity with the more or less 
monopolar network fields and lead to stronger shifts at the net- 
work b oundary. The r esult fo r our bright part data is similar to 
that of lAiouaz et al.l d2005bl) and the similar mechanism may 
also be responsible for it. But in the dark part, where coronal 
funnels were found dTian et al.l2008al) . the strongest blue shift of 
Ne vm appears to be in t he center part of the fu nnel. By using dif- 
ferent heating functions. lAiouaz et al.l d2005al) developed a two- 
dimensional MHD model. For one of the heating functions, the 
maximum blue-shift (outflow) was found not to be in the very 
center of the funnel but in the vicinity of the center, whereas, 
for the other heating function, the maximum was aligned well 
with the center of the funnel. Considering this model, the dif- 
ferent locations at which Ne vm reaches its maximum in the 
dark and bright part of our data set might be explained as the re- 
sult of different heating mechanisms dominating in the two parts. 
However, th e significant uncertain ties in both our measurements 
and those of lAiouaz et al.l d2005b ) imply that the above discus- 
sio n is not very co nvincing. 

lAiouaz! d2008l) claimed that the maximum inflow (redshift) 
at transition region temperatures is statistically toward the center 
of the network in the quiet Sun and toward the boundary of the 
network in the coronal hole, although large uncertainties were 
present in their plots. In our case, the redshift of the transition 
region line C iv increases steadily from internetwork to network 
center in both the bright a nd dark parts, which is similar to the 
result for the quiet Sun in lAiouaa d2008l) . The increase is more 
pronounced in the bright than the dark part. 

The Doppler shift of the chromospheric line Si n line remains 
at the zero level almost everywhere in and outside the network, 
indicating that the structures of chromospheric loops are similar 
everywhere in the lower solar atmosphere and there is a balance 
between upflow and downflow in small chromospheric loops. 

At the end of this section, we remind the reader of the prin- 
ciple defect of this dispersion plot method. As stated above, the 
distance to the network center is represented by, and the bins 
in Figs. [3] and [4] are divided according to, the filtered intensity 
of Si ii. A comparable intensity contrast between the network 
and internetwork for different supergranular cells is required of 
a successful statistical study of the radiance and Doppler shift 
distributions across the network. Although we applied the three- 
step filtering method three times to reduce local brightenings, 
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Fig. 2. Histogram of the filtered inten- 
sity of Si ii in the bright part (left) and 
the dark part (right). 
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Fig. 3. Dispersion plots for the bright part. The horizontal axis is the filtered Si n intensity {log(Ifat.si nl < Ifiit.si n >)), which can 
be considered to represent the distance from network center. A higher log(Im.si nl < Ifm.si u >) represents a position closer to 
the network center. The dashed line represents the outer edge of network boundary and the dotted line marks boundary of network 
center. 



there are still some relatively brighter and darker networks in 
Fig-Ill In Fig.Q] the green and yellow circles correspond to posi- 
tions of the data points in bins representing the network cen- 
ter in Figs. [3] and |4j respectively. It is clear that some data 
points at network boundaries were classified into the network 
center bins. The results in the network center bins were a mix- 
ture of contributions from some network centers and boundaries. 
In contrast, there are some networks that are insufficiently bright 
and have a lower filtered intensity of Si n. The intensities and 
Doppler shifts in the centers of these networks were not consid- 
ered as contributions to those of the network center bins. This 
mixture effect is difficult to remove. Applying the three-step fil- 
tering method repeatedly can reduce the intensity contrast be- 
tween the network and internetwork for different supergranular 
cells. It can also severely destroy the original network pattern. 
Our choice of three iterations should therefore be appropriate. 
However, the mixture effect is still existing and should not be 
neglec ted in our case. The a bov e discus s ion m ight also be the 
case in lAiouaz etall d2005bl) and lAiouazI d2008l) . 



5. Height variations of the extrapolated magnetic 
field 

In our previous works jTu et al.l2005btlHe et alj2007HTian et all 
l2008al) . we extrapolated the magnetic field from the photo- 
spheric magnetogram based on a force-free model for this data 
set. Here we used this previous result but studied the height vari- 
ations in the extrapolated magnetic field at different locations 
relative to the network center. 

For each of the two parts, we first selected three sets of grid 
points, which correspond to three bins (one bin at the network 
center, the other two bins at the network boundary and internet- 
work) of the filtered Si n intensity, at Mm in our calculation 
box. We then calculated the average values of the total magnetic 
field magnitude B and the ratio |J5 Z /J5|, where B z is the vertical 
component of B, in each set. After that, we selected three new 
sets of grid points which have the same x(east-west direction) 
and y(south-north direction) coordinates as those of the original 
sets of grid points at different heights in our calculation box. We 
then calculated the average values of B and \B Z /B\ in each new 
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Fig. 4. Dispersion plots for the dark part. The horizontal axis is the filtered Si n intensity (log(Ifu,.si nl < 1/ut.Si u >)), which can be 
considered to represent the distance from the network center. A larger log(Ifut.si nl < I/ut.Si n >) represents a position closer to the 
network center. The dashed line represents the outer edge of the network boundary and the dotted line marks the boundary of the 
network center. 
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Fig. 5. The height variations in B (left) 
and \B-JB\ (right) for the bright part. 
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set at differen heights. The height variations in B and \B Z /B\ in 
each set are presented in Fig. [5]for the bright part and Fig. |6]for 
the dark part. 

The left panels of the two figures clearly show that the mag- 
netic field is far stronger in the network than the network bound- 
ary and internetwork at Mm, which is a reflection of the well- 
known fact that magnetic flux tends to cluster in the network. We 
can see that the magnetic field strength decreases with height al- 
most exponentially at each of the three selected locations. The 
difference between the network and internetwork field also de- 
creases with height and reveals a more homogeneous magnetic 
structure in the upper part of the solar atmosphere. 

\B Z /B\ is an indicator of the inclination of a magnetic field 
line with respect to the horizontal direction. The right panels 
of the two figures reveal that \B Z /B\ is much larger inside the 
network than at the network boundary and internetwork region 
at lower heights above the photosphere. This is easy to under- 
stand because there are more open field lines and funnel struc- 
tures originating in the network than the internetwork. Above 
10 Mm, the difference in \B~/B\ between the three locations ap- 
pears, however, to decrease with height. In the bright part, they 
are equal to each other at about 30 Mm, where the value of \B Z /B\ 
is about 0.5. This result suggests a more homogeneous magnetic 
structure dominated by large loops with all possible i nclinations 
in the upper part of the solar atmosphere (Peter 2001). We found 
that the values of \B Z /B\ at 30 Mm are far larger (about 0.7-0.8) 
than 0.5 in the dark part, which is co nsistent with our previous 
result that there are open field lines dHe et al.ll2007t iTian et all 
l2008al) in this part. We note that the value of \B Z /B\ at 30 Mm 
above the network center is larger than values at the other two 
selected locations. This result may be explained that open field 
lines tend to stretch up vertically from the network centers. 

6. Average-profile study 

The second method used in our study of the distributions of radi- 
ance and Doppler shift across the network is an average-profile 
study. This method aims to study the average trend of radiance 
and Doppler shift from internetwork region to network by select- 
ing some individual cases. 

We chose 74 parallelogram regions, each of which stretches 
from one network cell to a neighboring cell across the bridging 
network; in our case they are close to the places where more 
than two network lanes intersect. In Fig. [7] the selected cases are 
outlined in black in the filtered intensity i mage of Si n. In both 
panels, the red network lines are taken from lHassler et al.l d 19991) 
and can be used as a good approximate of the network pattern. 

For each case, we first averaged the value of log(Ifm.si nl < 
Ifiit.si n >) along the short side of the parallelogram region and 
obtained the distribution of log{Ifu t snil < Ifm.siii >) along 
the long side. We found the position of the maximum value of 
log(Ifnt.si nl < 1/at.Si ii >) in the parallelogram region and de- 
fined this to be the network center, we then defined the two lo- 
cations where log(Ifut.si nl < I/m.si n >) declines to 92% of its 
maximum value to be the network boundaries (still in the paral- 
lelogram region) on both sides of the network center. By using 
this definition we can calculate the average size of the 74 net- 
works. We obtained the va lue of 9.88 ± 0.29 arcse c, which is 
consistent with the result of iPatsourakos et aTl (1 19991) . 

We then defined the network center to be the position 0, and 
the two network boundaries to be position 1. We extended the 
range 0-1 twice on both sides. The new boundaries were defined 
as 2. We divided the range 0-2 into 10 bins according to the posi- 
tion on each side, that is 0-0.2, 0.2-0.4, 0.4-0.6, 0.6-0.8, 0.8-1.0, 



Filtered intensity image of Si II (arbitrary unit) 




100 200 300 400 

X (second of arc) 

Doppler shift of Ne VIII (km/s) 




X (second of ore) 



Fig. 7. Upper: the filtered intensity image of Si n. The selected 
cases are outlined in black. Bottom: the Dopplergram of Ne vm. 
The large-blueshift cases are outlin ed in black. In both panels, 
the red network lines are taken fromiHassler et al] (I1999I) . 



1.0-1.2, 1.2-1.4, 1.4-1.6, 1.6-1.8, 1.8-2, which we referred to as 
the normalized distance from the network center. Then the two 
bins on different sides of the network center with the same nor- 
malized distance from the network center are combined into one 
bin. Therefore, there were only 10 bins for each case and we 
calculated the average values of intensity and Doppler shift for 
every line in each bin. 

We applied the method mentioned above to each case and 
obtained the profiles of the intensity and Doppler shift of every 
line across the network. We then normalized the intensity and 
Doppler shift in the form of A A ^' , where A is the value of line 
parameter(intensity or Doppler shift), and A m ,„ and A max are the 
minimum and maximum value of the parameter in the normal- 
ized distance range 0-2. We named this quantity the normalized 
difference of intensity or shift. For the Doppler shift of Ne vm, 
A is used as its minus value in order to study the blue shift of 
this line. We then obtained the normalized profiles of the line 
parameters across the network for each case. For all 74 cases, 
we extracted the 74 normalized values of each parameter in the 
bin at the same normalized distance, and calculated the mean 
value and the associated standard deviation. Then we obtained 
an average normalized profile of each line parameter, which is 
shown in Fig. [8] 

General trends of the radiance and Doppler shift distribu- 
tions across the network are revealed in Fig. [8] We found that the 
emissions of the three lines all declines from the network center 
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Fig. 8. Average profiles of radiance and Doppler shift for the selected cases. 



to internetwork region. The binned intensities of Si n and C iv 
decrease slowly from network center to network boundary and 
more significantly in the internetwork region. They have a corre- 
lation coefficient of 0.99. The correlation coefficient between the 
binned intensities of Si n and Ne vra is 0.95, indicating that the 
network structure is still prominent in the upper transition region 
of the quiet Sun. 

Although with large uncertainty, the general trend of the 
Ne vni blue shift decreases from the network center to the in- 
ternetwork region. The redshift of C iv is higher in the network 
and continues to decrease beyond network boundary. As for Si n, 
its Doppler shift appears to be slightly higher in the network than 
the internetwork, which may be due to the stronger magnetic ac- 
tivity in the network. These Doppler-shift curves can provide a 
measure of only the average trend and we are unable to draw firm 
conclusions because of the significant size of the uncertainties. 

7. Result of large-blueshift cases 

It has been found that large blue patches seen in the doppler- 
gram of Ne vm are associated with the chromospheric network 
dHassler et al. 1 I1999I) . Since a large blue shi ft of Ne vm is con- 
sidered to be a signature of the solar wind dHassler et alj 1999; 
Tu et al.ll2005ah or mas s supply to quiet coronal loops dHe et al.l 
20071: Tian et al.ll2008al) . it would be interesting to study the dis- 
tribution of the Ne vm blue shift in those network regions around 
where large blue shift of Ne vm is found. 

We selected cases for which the maximum Ne vm blue shift 
is higher than 5 km/s in the normalized distance range 0-2, from 
the 74 cases described in the previous section. We show 49 cases 
in the bottom panel of Fig. [7] for which our criterion is met. The 
mean size of the network is 9.05 + 0.53 arcsec for these large- 
blueshift cases. It is slightly smaller than the average network 
size of all cases. We then constructed in a similar way an aver- 
age normalized profile of each line parameters. The results are 
presented in Fig. [9] 

The contrast between the intensity in and outside the net- 
work is significant for Si n and C iv, which is natural and should 



be the case for almost every supergranular cell. However, the in- 
tensity profile for Ne vm is rather flat, which suggests that the 
magnetic fu nnels, whether being the footpoin ts of large loop s 
(lPeterll200lb or associated with the solar wind dTu et al.ll2005al) . 
expand and stretch into the region above the internetwork at the 
heights where considerable outflows of Ne vm are found. 

As for the Doppler shift, we found that the highest blue shift 
of Ne vm does not coincide with the highest redshift of C iv. The 
blue shift of Ne vm decreases steadily from the network center 
to the internetwork region, while the redshift of C iv reaches its 
maximum close to the network boundary (the normalized dis- 
tance 0.7-1.1). Thi s result appears to b e co nsistent with the s ce- 
nario proposed by lAxford et all (1 19991) and iTu et al.l (l2005cl) . in 
which reconnection occurs between the open funnels and side 
loops. The outflows produced by reconnections around a funnel 
tend to converge towards the center of the funnel. In contrast, 
the hot plasma trapped in low-lying loops are pulled down when 
they cool, and the downflows are stronger at the boundary of the 
network where side loops are accumulated. The Doppler shift 
of Si ii appears to be slightly higher in the network than in the 
internetwork region. We propose that it is due to the stronger 
magnetic activity in the network, although the structures of chro- 
mospheric loops are similar in and outside the network. 



The detailed distributions of radiance and Doppler shift 
across the network at different layers (chromosphere, transition 
region, and corona) above the photosphe re have not been stud- 
ied systematically . Our work with that of lAiouaz et al.l (2005b) 
an d lAiouazI (120081) may be a first step in analyzing these distribu- 
tions. The 1-D modeling of the acceleration of the solar wind has 
been investigated intensively. With the improvement of current 
observations, it is increasingly important to develop 2-D and 3- 
D models to study the origin and acceleration of the solar wind. 
We believe that the above observation results are helpful for 2-D 
modeling of network emission and outflows in coronal funnels. 
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Fig. 9. Average profiles of radiance and Doppler shift for the large-blueshift cases. 



8. Summary 

We have completed a statistical analysis of the distributions of 
both the radiance and Doppler shift across the network at differ- 
ent layers (chromosphere, transition region, and lower corona) 
above t he photosphere. We ap plied the filtering technique devel- 
oped bv lAiouaz et al.l (l2005bl) to an intensity image of the chro- 
mospheric line Si n obtained by SUMER/SOHO in a quiet Sun 
region. The filtered intensity image of Si n was used to reveal the 
network. Two different methods, namely the method of disper- 
sion plot and the average-profile study, were used to investigate 
the distributions of radiance and Doppler shift of Si n, C iv, and 
Ne vm across the network. 

We divided the entire quiet Sun region into a bright and a 
dark part, and produced dispersion plots for each part. The in- 
tensities of Si n and C iv show an obvious increasing trend from 
network boundary to network center in both parts. However, the 
intensity of Ne vm shows different trends, namely increasing in 
the bright part and decreasing in the dark part. This result might 
suggest a stronger expansion of magnetic structures in the dark 
part than in the bright part. In each of the two parts, the Doppler 
shift of C iv increases almost steadily from internetwork to net- 
work center. The average Doppler shift of the chromospheric line 
Si ii remains at the zero level in and outside network. The distri- 
bution of Ne vm Doppler shift is not clear due to the large uncer- 
tainties. This method largely depends on the intensity contrasts 
between the network and internetwork for different supergranu- 
lar cells. The mixture effect influences the statistical results and 
could not be neglected because it is difficult to obtain a compa- 
rable contrast for different supergranular cells in real conditions. 

The difference between height variations of the extrapolated 
magnetic field above the internetwork, network boundary, and 
network center was also investigated; this difference implies 
a more homogenious magnetic structure at greater height and 
clearly reflects the different magnetic structures in the two parts. 

By selecting individual cases and investigating the average 
profiles of radiance and Doppler shift distributions across net- 
work, we found that the intensities of the three lines all decline 



from the network center to internetwork region. The correlation 
coefficient between the binned intensities of Si n and Ne vm is 
0.95, indicating that the network structure is still prominent at 
the layer where Ne vm forms in the quiet Sun. Although with 
large uncertainties, the average Ne vm blue shift decreases from 
the network center to internetwork region. The redshift of C rv is 
higher in the network and decreases beyond the network bound- 
ary. 

Since the high blue shift of Ne vm could be a signature of 
solar wind or mass supply to coronal loops, we also studied the 
average profiles of radiance and Doppler shift across network 
for the large-blueshift cases. The intensity distribution of Ne vm 
does not correlate with those of Si n and C iv. We found that the 
blue shift of Ne vm decreases steadily from the network center 
to the internetwork region, and the large blue shift of Ne vm does 
not coincide with large red shift of C iv. 

Our observation results are worthy of consideration in terms 
of 2-D modeling of network emissions and outflows in coronal 
funnels. 
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